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Abstract 

Parametric enhancement of the appearance probability of the neutrino os- 
cillation under the inhomogeneous matter is studied. Fourier expansion of 
the matter density profile leads to a simple resonance condition and man- 
ifests that each Fourier mode modifies the energy spectrum of oscillation 
probability at around the corresponding energy; below the MSW resonance 
energy, a large-scale variation modifies the spectrum in high energies while 
a small-scale one does in low energies. In contrast to the simple parametric 
resonance, the enhancement of the oscillation probability is itself an slow os- 
cillation as demonstrated by a numerical analysis with a single Fourier mode 
of the matter density. We derive an analytic solution to the evolution equa- 
tion on the resonance energy, including the expression of frequency of the 
slow oscillation. 
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Experiments over the last decade have established the oscillation among 
the neutrino flavors and narrowed the allowed value of oscillation parameters. 
Experiments of the next generation aiming for precision measurements are 
actively discussed, and many of them involve a neutrino oscillation under a 
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very long baseline penetrating through the Earth. The matter density on 
such a path varies spatially, reflecting the interior structure of the Earth. 
The neutrino oscillation under the varying density of matter has been widely 
studied in works including Refs. [l], 0, HJH, and the oscillation probability 
is found to be significantly enhanced in some cases. This enhancement can 
be understood in a context of the parametric resonance 0, 0] . These studies 
usually model the density profile by the step function segmented into three 
pieces, among which the middle piece has a larger value than the both ends, 
and reproduces the qualitative feature of the numerical results 0]. In this 
paper, we present an alternative formulation based on a Fourier expansion 
of the matter density profile, which generalizes the sinusoidal matter profile 
studied in one of the earliest studies [2J. This approach organizes the effects 
of the matter profile on the neutrino oscillation by their frequencies, in terms 
of which the condition for the parametric resonance to take place is given. 

We will analyze the neutrino oscillation in the two-flavor framework rather 
than the standard three-flavor one to make the discussion simple and clear. 
We briefly explain that this simplification does not wipe the effect of the 
matter density profile off the oscillation probability. Figure [T] shows the 
v c appearance probability for the baseline length L = 12 000 km. 
The curves (i) and (ii) are calculated in the three-flavor framework. The 
curve (i) assumes the matter density p(x) that varies according to the Pre- 
liminary Reference Earth Model (PREM) [5|, while (ii) assumes the con- 
stant density of 7.58 g/cm 3 , which is the average value of the PREM profile. 
Other values of parameters taken here are Sm^ = 7.9 x 10~ 5 eV 2 , = 
2.5 x 10- 3 eV 2 , sin 2 26 12 = 0.84, sin 2 2023 = 1.0, sin 2 26 13 = 0.05, and 5 C p = 0, 
in the standard notations of Ref. We find qualitative difference between 
(i) and (ii), showing a major effect of the inhomogeneity of the matter den- 
sity. The curves (iii) and (iv) are obtained from the two-flavor calculation 
renormalized by sin 2 #23- The curve (iii) assumes the cosine density profile 
of the matter as [7.58 — 4.32 cos(27rx/L)] g/cm 3 , where the two values are 
taken from the first two Fourier coefficients of the PREM profile. The curve 
(iv) assumes the constant density of 7.58 g/cm 3 . The mass-square difference 
and the mixing angle are taken in common as Sm 2 = 2.5 x 10~ 3 eV 2 and 
sin 2 29 = 0.05. These two curves again show the significant qualitative dif- 
ference from each other. Meanwhile, we find the curves (i) and (iii) are quite 
similar except in the low energies, and so are (ii) and (iv). This suggests the 
following two points: first, the two-generation analysis is sufficient to analyze 
the essence of the matter profile effect; secondly, the Fourier analysis of the 
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Figure 1: Oscillation probability P(v^ — > v c ) for L — 12 000 km. The probabilities (i) and 
(ii) are evaluated by three-generation analysis and (iii) and (iv) by two-generation analysis. 
The matter density profile used for (i) is the PREM profile [5J], for (ii) and (iv) is the 
constant profile with 7.58 g/cm 3 , for (iii) is the profile with [7.58 — 4.32 cos(27ra/L)] g/cm 3 . 
The values of other parameters taken are described in the text. 
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matter density profile is effective in analyzing its effect. 

We study the parametric resonance of the two-flavor neutrino oscillation 
from Ufj, to u e . The evolution equation in matter is given by 
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where E is the energy of neutrinos and a(x) = 2y/2EGFNAY e p(x)/[1.0 g/mo\] 
accounts for the matter effect. Here Gp, and Y e are the Fermi constant, 
the Avogadro constant, and the proton-to-nucleon ratio, respectively. We 
reduce Eq. flTJ to a differential equation of v e (£) m a dimensionless form as 



+ iA m (£K(£) 
-A 2 A cos 29 
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where £ = x/L with baseline length L, A m (£) = a(£)L/2E, and A = 
5m 2 L/2E. We eliminate the first- derivative term using 



ds iA r 



z(0 = z/ c (£)exp 

and then expand the matter density into a Fourier series as 
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where po is the average density [lj. Accordingly we expand 
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to obtain 
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which we read as an equation of a harmonic oscillator, where 



uo 2 = - (A m0 - Acos2#) 2 + A 2 sin 2 20 + D (6b) 

with D = J2^Li gi ves the natural frequency of the system, while 

M n(0 = cos 2n7r£ — if3 n sin 2t7,7t^ + 7„ cos 4mr£, (6c) 

vni(0 = ™ ml [cos2(n+l)irt + cos2(n-l)ir£\, (6d) 

with 

a n = A m „(A m0 -Acos2#)/2, (7a) 

/3 n = n7rA mn , (7b) 

In = AL/8 (7c) 

give oscillating potentials. The initial conditions v e (0) = and ^(0) = 1 
become z(0) = and z'(0) = -i(A/2) sin 29. 

We recall that a harmonic oscillator with an oscillating potential is led to 
a parametric resonance when the frequency of the oscillating potential is, in 
the simplest case, twice the natural frequency of the system [7J. Applying this 
resonance condition to our system, we find that the n-th Fourier mode of the 
matter profile, whose frequency is 2mr, gives rise to the parametric resonance 
of the neutrino oscillation when 2mx = 2uj. This resonance condition reads 
E = E { n ±] with 

E'f> = . Sm2L/2 , (8) 



A m o cos 

29 ± ^/An 2 7T 2 - A^ sin 29 — D 



where E n + ^ lies below the MSW resonance while E n ' lies above. The energy 
spectrum of the appearance probability at this energy hits the local minimum, 
which we call n-th dip, as is evident from the initial condition u e (0) = and 
the condition uj = nn. To summarize, each Fourier mode of the matter 
density is expected to resonate the appearance probability of neutrinos at 
around the n-th dip of the energy spectrum. 

We elucidate the effect of the n-th Fourier mode by considering a simpli- 
fied model where the matter density profile consists only of the n-th mode 
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around the average density. The matter effect in this case is described by 
A m (£) = A m0 + A mn cos2n7r£, and Eq. fl§]) is reduced to 



which is free from the v n i terms. This is Hill's equation, which is a general- 
ization of Mathieu's equation that appears in typical parametric resonance 
problems . To verify the discussion in the previous paragraph, we numeri- 
cally solve Eq. (Q and show in Fig. [2] the oscillation probability Piv^ —>■ u e ) 
for the various values of p\ in (a), p 2 in (b), and p 3 in (c). The average 
density is fixed to po = 7.58 g/cm 3 and other parameters are taken same as 
for (iii) and (iv) in Fig. [U The graphs indicate that the first, the second, 
and the third Fourier mode of the matter density profile indeed modify the 
appearance probability around the first, the second, and the third dip, at 
which ojq equals to n, 2ir, and 3n, respectively. 

We further confirm explicitly the presence of the resonant behavior. For 
simplicity, we take the energy of the neutrino at the resonance energy E n + ^ 
below the MSW resonance, while the following analysis is equally applicable 
to the case of E = En . We consider the propagation of neutrinos through 
a repetitive density profile and calculate P{y^ — > v e ) as a function of £ 
up to £ > 1 so that the resonant behavior becomes evident. The result is 
shown in Fig. [3], where (a) p\ = — 0.3g/cm 3 , (b) p\ = — 1.0g/cm 3 , and (c) 
P2 = —0.3 g/cm 3 . The energy of neutrino is accordingly fixed at e[ + ^ in (a) 
and (b) and E% in (c), while other parameters are taken common and same 
as in Fig. [2J The dotted curves, evaluated by the analytic formula of Eq. ffTS]) . 
will be discussed later. We find that the repetitive matter profile indeed raises 
the oscillation probability as expected from the parametric enhancement. 
Unlike the resonance in the classical dynamics, however, the unitarity bound 
prevents the probability from the divergence. The enhanced probability turns 
to decrease after reaching its peak, and thereby lies on a slow oscillation so 
that the probability consists of the fast and the slow oscillation. 

We derive an analytic formula of this parametric-enhanced oscillation. 
We still keep E = E n + ^ here so that ujq = nir and A = 5m 2 L/2E^\ and 
solve Eq. ([9]) as 



z "(£.) + ( w o + a n cos 2n7r£ — \(3 n sin 2n7r£ 

+ j n cos4n7r£)z(£) = 



(9a) 




(9b) 



z(£) = A(£) cosniit; + -8(0 sinn7r£ , 
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Figure 2: Oscillation probability Piv^ — > f c ) for p„ = 0.0 g/cm 3 , —0.5 g/cm 3 , —1.0 g/cm 3 , 
and —3.0 g/cm 3 with (a) n = 1, (b) n = 2 and (c) n = 3. The average density is fixed to 
po = 7.58 g/cm 3 and other parameters are taken same as for (hi) and (iv) in Fig. [TJ 
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Figure 3: Oscillation probability P{y^ — ► v c ) beyond £ = 1. The solid curve is given by 
the numerical calculation and the dotted one by using the formula of Eq. (fT5|) . The matter 
profile consists of p\ — —0.3 g/cm 3 in (a), pi = — f.0g/cm 3 in (b), and p 2 = —0.3 g/cm 3 

in (c). Neutrino energy is fixed at e[ + ' in (a) and (b) and E^ in (c). The values of other 
parameters are same as in Fig. [2j 
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where A(£) and B(£) depend weakly on £. We then have 

cos mrl; 



+ 



sin mr^ 



:n) 



+ (faster-oscillating terms) = . 



We neglect the faster-oscillating terms and assume A(£),B(£) oc e lf ^ in ac- 
cordance with Floquet's theorem [j| to obtain 
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Hence Q is given by 
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which is real since ujq 



/3 n /A mn and /3 2 
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(13) 



AL(2A 2 sin 2 20 + 



A 2 nn )/8 > 0. Thus A(£) and oscillate with the frequency Q, and give 
z(£) the slow oscillation mentioned above. The oscillation probability can be 
enhanced by this slow oscillation while it does not diverge due to its unitarity, 
in contrast to usual parametric resonances. When \J f3 2 — a 2 / (Amr) < 1, we 
obtain 

rm _ a 2 a 



A 2 sin 2 20 + \^ mn 



Ann 8rm V ' 2~ mn " ^ 

Taking the lowest order in A mn of Eq. (fill) , we reproduce f2 n in Ref. [2| with 
2c<} = nu and n = 1. 

The solution of Eq. with the initial conditions u c (0) = and ^(0) = 1 
is expressed in terms of the oscillations with the frequencies uq = nn and Q 

as 



iA sin 26 
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We calculate |^(0| 2 from Eq. (fT5|) and plot in Fig. E] by the dotted curves. 
There we see that the frequency of the slow oscillation becomes larger for 
smaller n and larger \p n \, as we can read from Eq. (TH1) . The frequencies of 
the oscillation, both of the fast and slow ones, agree well with the numerical 
results. The amplitude of the slow oscillation agrees well with the numerical 
result in (a) and (c), while the disagreement between the two grows as the 
frequency Q becomes large in (b). 

The slow oscillation is extracted by averaging out the fast oscillation with 
the frequency rnr as 



A sin 29 and Q <C a n , f3 n <C 1 are satisfied. This result agrees with Eq. (2) of 
Ref. [2J in the lowest order in A nm with 2uo = nuj and n = 1. 

We studied the parametric resonance in neutrino oscillation under the 
inhomogeneous matter density. We expanded the matter density profile into a 
Fourier series and expressed the evolution equation of the neutrinos in a form 
of Hill's equation, which leads the oscillation probability to the parametric 
resonance. Noting that the resonance energy is controlled by the wavelength 
of the spatial variation of the matter density, we pointed out that the n- 
th Fourier mode of the matter density modifies the appearance probability 
through the resonance at the energy around the n-th dip. We verified the 
correspondence between the Fourier mode and the modification of the energy 
spectrum of the oscillation probability. 

We further confirmed the parametric enhancement of the oscillation prob- 
ability in a calculation where the matter density profile is periodically re- 
peated. Unlike a simple parametric resonance, the enhancement turned out 
to be a slow oscillation, which respects the unitarity. We developed an an- 
alytic formula of the oscillation probability under the parametric enhance- 
ment. It well reproduced the frequencies of the slow- and the fast-oscillation 
of the numerical results. The agreement on the amplitude of the slow oscil- 
lation diminishes as the frequency of the slow oscillation gets larger. 

The Fourier analysis of the parametric resonance cleared up the corre- 
spondence between the energy of the neutrinos and the scale of the spatial 
variation of the matter density. Below the energy of the MSW resonance, 




(16) 



This equation reduces 




when the relations A 
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a large-scale variation modifies the oscillation spectra in high energies while 
a small-scale one does in low energies. The effect of the density variation 
becomes large when neutrinos propagate through the core of the Earth. This 
effect should be kept under control in analyzing the oscillation experiments 
with very long baselines, including the next-generation accelerator experi- 
ments and the atmospheric neutrino ex per iments. Conversely, it can serve 
as a probe for tomography of the Earth [lOj once the resonance effect is well 
understood. Systematic Fourier analyses on the parametric resonance will 
give a perspective to these experiments. 
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